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We study the influence of filament elasticity on the motion of collective molecular motors. It is 
found that for a backbone flexibility exceeding a characteristic value (motor stiffness divided through 
the mean displacement between attached motors), the ability of motors to produce force reduces as 
compared to rigidly coupled motors, while the maximum velocity remains unchanged. The force- 
velocity-relation in two different analytic approximations is calculated and compared with Monte- 
Carlo simulations. Finally, we extend our model by introducing motors with a strain-dependent 
detachment rate. A remarkable crossover from the nearly hyperbolic shape of the Hill curve for 
stiff backbones to a linear force-velocity relation for very elastic backbones is found. With realistic 
model parameters we show that the backbone flexibility plays no role under physiological conditions 
in muscles, but it should be observable in certain in vitro assays. 



o 



I 

S3 
O 

o 



(N 
> 

(N 



ON 



T3 

S3 
O 

o 



I. INTRODUCTION 

Molecular motors play a key role in a variety of bi- 
ological processes like muscle contraction, intracellular 
transport, cell locomotion, flagellar rotation etc. M. De- 
spite structural similarities, motors can be classified into 
two groups according to their function. Processive mo- 
tors 0, also called "porters" ||, consist of a single 
molecule which can move over long distances along its 
molecular track without dissociating from it. The most 
common processive motor is kinesin interacting with mi- 
crotubules. Nonprocessive motors, also called "rowers", 
can only generate macroscopic motion when operating in 
large groups. Muscular myosin, interacting with actin, 
belongs to this class of motors. Here we focus on non- 
processive motors. 

For many decades exclusively data from physiologi- 
cal measurements on muscles Q provided experimen- 
tal information for modeling collective molecular motors 
[^|,^). In recent years, a variety of in vitro experimen- 
tal techniques allowed the observation of single motor 
proteins. These experiments include gliding assays 
optical tweezers |),|l0| and micromechanical force mea- 
surements flllfl . They allowed for a new insight into the 
basic principles underlying the operation of motors. Not 
only new theoretical models for single-molecule motors 
Jl2| |l5| were inspired by these experiments, but also new 
models for cooperative motors P, ^6|Jl7| . Except for the 
work by Csahok et al. [Q, which discusses the transport 
of elastically coupled particles driven by colored noise, 
all these models deal with motors placed on a rigid back- 
bone, interacting with a rigid track. 

The assumption of stiff filaments seems to be appropri- 
ate under physiological conditions in muscles, since the 
measured extensions of few nanometers jl9] ^lj are suffi- 
ciently small compared to the myosin step-size, which is 
about 10 nm p^ j. However, it certainly can become in- 
valid in gliding assays of the type discussed in Refs. J7||23[ 



if an actin filament glides over widely separated linearly 
placed motors. If the spacing between motors is large 
enough, the elasticity of the backbone or track section 
between two motors can become comparable to the elas- 
ticity of a single motor head. Experiments with myosin 
molecules bound to an elastic background are conceivable 
as well. 



It is the purpose of this paper to investigate qualita- 
tively and quantitatively the influence of filament elas- 
ticity on the operation of myosin-like motors. Our major 
quantity of interest is the force- velocity-relation (filament 
velocity as a function of the external load) . We want to 
identify the universal effects of filament elasticity and at 
the same time keep the model as close as possible to ex- 
periments. The elasticity modeled by linear harmonic 
springs may either originate from the flexible backbone 
or from the flexible track (Fig. |l|). As long as we are 
dealing with small relative elongations, both sources of 
flexibility are equivalent and our model should apply to 
both cases. For specificity we use a formulation with an 
clastic backbone and a stiff track. 





FIG. 1. Two possible sources of elastic coupling of col- 
lective molecular motors: elastic backbone (upper figure) 
and elastic track (lower figure). Except for very soft back- 
bones/tracks, the description of both models is equivalent. 
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Another important component of the model is the 
modeling of individual motor heads. Most models de- 
scribe the heads by several states with different confor- 
mations and transition rates between them. Such models 
are based on early ideas by A.F. Huxley || and have later 
been refined in order to explain more experimental de- 
tails. More recently simplified models using a two-state 
ratchet formalism have been developed in order to con- 
centrate on generic features of motion generation Jl7| . 
Another class of models describing Brownian particles 
in ratchet potentials subject to colored noise, however, 
seems to be only very distantly related to motor pro- 
teins p4[ ]. We decided to use a two-state crossbridgc 
model, similar to the model introduced by Leibler and 
Huse ||, but with just two long-living states. It includes 
the transitions between attached and detached state and 
an active power stroke, which have both been identified 
as basic elements of the myosin motor |||| an d also ob- 
served in vitro fTofl . Compared to other investigations 
P,p|JT7|] our model also contains a low number of free 
parameters, which makes it more suitable for a study 
of universal aspects. Yet, adding the strain-dependence 
of the detachment rate, the model we use is sufficient 
to describe the experimentally measured force-velocity- 
relation of actin/myosin in muscles ||J. A more detailed 
description of this two-state model can be found in Ref. 

HI- 

Due to the generic features of such two-state models we 
expect that the effects discussed here should as well ap- 
ply to other models which contain_the same basic mech- 
anisms of force generation, e.g. 
of the type discussed in Ref. flic- 




However, models 



which describe par- 
ticles in an asymmetric periodic potential subject to a 
temporally correlated noise, are based on a thoroughly 
different driving mechanism. Therefore, one expects and 
actually finds a variety of disparate effects including a 
strong influence of the coupling strength on the velocity 
even without external load. 

The outline of the paper is as follows. In Sec. [n] wc 
introduce the model ]26| ], describe its phenomenological 
properties, present the main results and discuss its im- 
plications to experiments. The full calculation for strain- 
independent detachment rates is shown in Sec. Ill and 
for strain-dependent rates in Sec. IV. 



II. DISCUSSION 

A. Description of the model 

We consider a one-dimensional model describing many 
motors which produce force between two filaments glid- 
ing past each other. The force is generated by a confor- 
mational change ("power stroke") in the molecular mo- 
tors fixed to the backbone, which takes place after they 
attach to the molecular track (Fig. ||). We assume that 
the motor proteins can be found in two states: attached 



to and detached from the track. This corresponds to 
taking into account only the two long-living states in the 
model used by Leibler and Huse [|| . The transitions be- 
tween these two states occur stochastically. We denote 
the mean life time of the attached state by t a and of the 
detached state by td- Each attached motor is described 
as a harmonic spring connecting its root at the back- 
bone (position y) and its head on the track (position x). 
The force this motor produces between the track and the 
backbone is given as k(x — y), with a spring constant k. 
Since the motor is in a forward leaning position before 
attachment (see Fig. ^), it attaches to the track at the 
point x n , which is the position y of the root of the motor 
before attachment, shifted by the displacement d ("power 
stroke" ) 

x n = y + d. (2.1) 

We neglect thermal fluctuations and the discreteness 
of binding sites on the track. This is motivated by our 
recent results for rigid backbones |25|] , where we find that 
discrete binding sites (with a spacing of 5.5 nm) as well 
as thermal fluctuations have only a minor effect on the 
resulting force- velocity-relation. 

(i-i) (i) (i+i) 

(a-l)(a)(a+l) 
Za L 

Backbone 



Track 



FIG. 2. Definition of the model. Motors are fixed on the 
elastic backbone at uniform spacing L and attach to the stiff 
track. z a denotes the position of a-th motor on the unstrained 
backbone. Due to the conformational change, the head of each 
motor attaches at the distance d from its root. 

While assuming a stiff track, we model the backbone as 
a linear spring with compliance 7 _1 per unit length. But 
note that this is merely a convention. Our results apply 
equally well to the reverse situation, where the molecular 
motors are fixed on a (rigid) cover slip and the elasticity 
is due to the molecular track transported by them. We 
consider N motors placed on the backbone at uniform 
spacing L, so that the total backbone length is NL. Note 
that the assumption of uniform spacing is made solely for 
simplicity. Any other distribution which is homogeneous 
on length scales L would lead to t he s ame exponential 
distribution of gap widths (see Eq. 2.3 below). The po- 
sition of the a-th motor on the unstrained backbone will 
be referred to as z a . In the following, the actual posi- 
tions of motor heads x a (t) and of motor roots y a (t) are 
measured relative to z a (Fig. ||). 

Instead of using the quantities N and L it will prove 
helpful to use the mean number of attached motors 
N = Nt a /(t a + td) and the mean spacing between two 
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attached motors L = L(t a + td)/t a . Since for now we are 
dealing with strain-independent reaction rates the distri- 
bution of attached and detached states does not depend 
on the motion of their positions. The probability of find- 
ing a gap with a detached motors between two attached 
ones is given by the geometric distribution 



t a 



ta+tdj ta + td 



(2.2) 



As shown in Sec. Ill we find that in case of a time- 



independent external force it does not matter whether 
this force pulls on one end or homogeneously on the whole 
backbone with a density / = F/NL. A quite remarkable 
result of our analysis is that the force- velocity relation re- 
mains linear for flexible backbones and that the zero load 
velocity d/t a does not depend on the backbone elastic- 
ity. The force-velocity relations for a stiff and an elastic 
backbone are shown in Fig. ^. 



FIG. 3. The moving backbone (solid line) and its initial 
position at time t = (dashed line). y a denotes the root po- 
sition of the a-th motor, relative to its initial value (z a ). x a 
denotes the head position of the a-th motor, also relative to 
the initial position of the a-th root, z a . 

In order to keep the model as lucid as possible we 
assume a small duty ratio [Q, meaning that a motor 
molecule spends most of its time in the detached state, 
t a *C td- Since we are dealing with nonprocessive motors 
("rowers" Q), this assumption is certainly valid. While 
keeping the mean number N of attached motors and their 
average spacing L constant, we consider the limit L — > 0. 
With this simplification the model becomes continuous. 
Also, the assumption about equidistantly placed motors 
on the backbone becomes superfluous in this limit. The 
distribution of motors is insignificant as long as it is suf- 
ficiently homogeneous on the length scale L. The at- 
tachment rate per length L (between the positions z and 
z + L on the backbone) obeys r a = L/Ltd = l/t a - The 
distribution of gap widths ( |2.2| ) takes the form of an ex- 
ponential distribution 



p(l) = \e-^ 



B. Results 



(2.3) 



In this subsection we summarize our main results for 
the analysis of the model described above; the details of 



the calculation are given later in Sec. [II 



As described by now the model contains seven indepen- 
dent parameters: t a , N, L, F, k, 7 and d. Upon measur- 
ing the force per motor, F/N, in units of the force during 
one power stroke, kd, and measuring the backbone elas- 
ticity per unit length, j/L, in units of the motor head 
elasticity, they may be reduced to two adimensional pa- 
rameters: F = F/Nkd and 7 = 7/fcL. Then the velocity 
in units of a single motor velocity v = vt a /d is given by 
a "scaling" function v — i](F,j). 



vt /d 




FIG. 4. Force-velocity-relation for the stiff (7 = 00) and 
elastic (7 = kL) backbone. 
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FIG. 5. The slope of the force-velocity relation (coefficient 
C2 in Eq. (3.16)), as a function of the relative backbone stiff- 
ness -y/kL: MC-simulations (stars), analytic approximation 
without correlations (lower curve), with correlations to the 
distances to nearest neighbors (upper curve). In both ana- 
lytic curves the numeric result for k(l/K e ) from Fig. was 
used. 



If the relative backbone stiffness 7 lies below 1, the 
slope of the force-velocity curve and consequently the 
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stall force differ significantly from those for a stiff back- 
bone (Fig. [sj) . While the stall force / st aii is proportional 
to the motor stiffness (/stall = kd/L) for stiff backbones, 
it becomes a function of backbone stiffness for very flex- 
ible backbones and is given as l^djvl?. Here v is a 
numerical constant which has the value v f .64 as ob- 
tained from the Monte-Carlo simulation. Using analyti- 
cal tools based on a Master-Equation approach with cor- 
relations between the position of a motor and the dis- 
tances to its neighbors, we obtain the value v w 1.50, 
which is in good agreement. For completeness Fig. |5| 
also shows the result of a Master-Equation without cor- 
relations, yielding v = 1 . 



Finally, we extend the model described in Sec. 11 A by 
introducing a strain-dependent detachment rate. This 
extension is inevitable for a quantitative comparison with 
experiments on the actin-myosin system. We already 
mentioned that (with or without backbone elasticity) 
strain-independent transition rates lead to a linear force- 
velocity-relation. However, since the very beginning of 
muscle research it has been known that the force- velocity- 
relation rather has a hyperbolic form, also called the Hill 
curve Q . ft has also been known that the energy libera- 
tion in a stretching muscle depends on velocity, which is 
also called the Fenn effect 27 1. The most natural way to 



reproduce these physiological measurements is the intro- 
duction of a strain-dependent detachment rate, meaning 
that the lifetime of the attached state t a = t a (xi — yi) 
is larger for those heads that have just gone through the 
power stroke and produce maximum force than for those 
which have already done their work and now pull back- 
wards. As a consequence the duty ratio becomes lower 
at higher velocities. This idea has already been used by 
A.F. Huxley ||. Although there is some direct exper- 
imental evidence for strain-dependent detachment rates 
JTot , the functional form of this dependence has not been 
measured yet. For simplicity we model this dependence 
as an exponential t a (£) = exp(a£), which suffices to fit 
the force- velocity-relation jl. However, we stress that 
this is only a first approximation and that other forms 
are possible as well. Further experimental information is 
highly desirable for a future more detailed modeling of 
molecular motors. 

Some other functional forms of the detachment rate 
(e.g. i Q (£) of £~ 2 ) lead to anomalous force-velocity- 
relations already with rigid backbones |25j ]. These can 
lead to oscillations, similar to those proposed by Jiilicher 
and Prost psf . We expect that flexible backbones can 
lead to additional phenomena like wave generation. 

Strain-dependent detachment rates enhance the diffi- 
culty of an analytic solution of our model enormously, 
since the distribution of attached and detached mo- 
tors depends on the distribution of head positions (xi). 
Therefore, we will mainly use Monte-Carlo simulations 
and restrict analytic arguments on limiting cases. The 
simulations show that two major analytic results of 
the strain-independent case carry over to the strain- 
dependent case: (i) If the backbone flexibility exceeds 



its characteristic value, the stall force decreases strongly, 
(ii) The backbone flexibility has only little influence on 
the zero-load velocity. 

For stiff backbones (j/kL ^S> I) the force- velocity- 
relation as measured by Hill B is reproduced perfectly 
using the value ad — 0.58. For decreasing backbone stiff- 
ness (j/kL < 1) the stall force decreases and the force- 
velocity-curve becomes increasingly linear. This can be 
understood as follows: The forces which lead to positive 
velocities become smaller and so does the strain on single 
motors. The small strain does not have any significant 
influence on the detachment rate any more and the re- 
sults obtained for a strain-independent detachment rate 
become exact. The crossover from the Hill curve to a lin- 
ear relation is shown in Fig. [| There ad = 0.5 is used, 
but the curves would look qualitatively similar for other 
positive values. 



vt/d 
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FIG. 6. Monte-Carlo results for force-velocity-relations 
with strain-dependent detachment rates (ad — 0.5) and dif- 
ferent backbone stiffnesses. Note the crossover from the Hill 
curve at stiff backbones to a nearly linear relation at soft 
backbones. 



C. Implications for experiments 



In order to apply our theory to experiments, we need 
the spring constants of the myosin heads and of the 
actin filaments. The elasticity of the attached myosin 
head was measured by Finer et al. |^] with the result 
k = 0.4pN/nm (measurements by Ishijima et al. p9| ] 
yield k — 0.28 pN/nm, which is in quite good agreement). 

In our model we assumed springs with flexibility pro- 
portional to their length, obeying fc spr i n g = j/L. Actu- 
ally actin is a semiflexible polymer. Its elastic behavior 
was subject of many theoretical and experimental stud- 
ies in last years [p0|-p4[ . There are essentially two contri- 
butions to the elasticity of actin: the longitudinal clas- 
tic modulus and the buckling of the polymer, induced 
through thermal fluctuations. 
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At very high loads, the stiffness is limited by the elastic 
modulus of actin and is proportional to L _1 . For a fila- 
ment of length 1 fim the stiffness is about 44pN/nm j34|. 
In this regime, the characteristic distance between at- 
tached motors can be estimated as L c h = j/k ~ 100 //m. 
In this case backbone flexibility is of no experimental rel- 
evance. 

At low tension, the buckling modes limit the stiffness of 
actin. In the linear response approximation, the stiffness 
of a polymer with length L obeys the law 



tpolymcr — QOksTipL 



(2.4) 



where l v denotes the persistence length. Although we 
assume springs obeying an L _1 -law, we can still use the 
L _4 -law to give an estimate for the characteristic dis- 
tance for which fc po i ymcr (L c h) ~ k. 

Recent measurements provide the value £ p — 7.4/im 
p5| . With these values we finally estimate the char- 
acteristic distance between attached motors L c h = 
(90/c B T£2/fc) 1 /4 ~ 500 nm. If the mean displacement be- 
tween attached motors is larger than L c h, we expect that 
the effects of backbone elasticity should be observable. A 
simulation with springs obeying the L~ 4 -law ( [2.4| ) is in 
preparation. 

In muscles 500 nm is about the length of a half sar- 
comere [p6P]. A rough estimate (300 myosin heads in 
one thick filament, 3 actin filaments per one thick fila- 
ment, 2.5 — 10% of heads in the attached state) leads 
to L values between 50 nm and 200 nm, significantly be- 
low the characteristic length L c ^. This implies that the 
elasticity of actin filaments does not influence the opera- 
tion of muscles. This result is not surprising - backbone 
elasticity always reduces the efficiency of motors and it 
would be hard to understand why muscles spoil their high 
efficiency in such a prodigal way. 



III. ANALYTICAL SOLUTION 

In this section the calculation leading to the force- 
velocity-relation for strain-independent detachment rates 
is given. This is done in several steps: first we show that 
the model behaves equivalently if a force acts on one 
end of the filament or homogeneously along the whole 
length. As second we calculate the effective compliance of 
a semi-infinite chain, which is an important input quan- 
tity for later use. Then we show the linearity of the 
force-velocity-relation, show that the zero-load velocity 
does not differ from its value for stiff backbones and fi- 
nally calculate its slope in two different approximations. 

A generic situation found in many experimental setups 
is that a force F, usually produced by an optical tweezer, 
acts on the rear end of the backbone. Another possibility 
is to produce the force by a viscous liquid, flowing along 
the backbone. Such a force acts more or less homoge- 
neously on the whole backbone. In both cases the main 



quantity of interest is the resulting mean backbone ve- 
locity depending on the load F. In the following we show 
that both situations are equivalent within the scope of a 
theoretical description. 



A. Point force 

From now on we use the index i, which runs over at- 
tached motors only, instead of the index a, running over 
all motors. Xi and yi denote the head and root positions 
of i-th attached motor relative to its initial position on 
the unstrained backbone {zi) at time t = 0. The stiffness 
of the backbone fragment between the motors i and i + 1 
equals 7/(zj+i — Zi). At the point where the i-th motor 
is fixed to the backbone, the sum of the all three forces 
(from the motor, from the left part of the backbone and 
from the right part of the backbone) must be zero: 



k(xi - yi) - 7 



Vi ~ Vi-i 



Zi 



+ 7 



K x i - Vi) +7 



Vt+i - Vi 

Zi+1 — Zi 
2/2 



yi 



o 



= F 



(3.1a) 
(3.1b) 



Z2 - Z\ 

The second equation describes the first attached motor 
and differs from the others since the backbone force act- 
ing from the left is replaced by the external force F. With 
given Xi and Zi this set of equations allows us to deter- 
mine the values of j/,. 

The detachment rate equals t~ l for each motor. The 
detachment of one motor is described by canceling its 
position in the set of x- and z-values. Afterwards, all the 
y- values are determined anew from Eq. ( |3.l| ). 

The process of attachment occurs at the rate N/t a 
and consists of choosing a random position z n between 
and NL, calculating the corresponding y(z n ) (the 
root position of the new motor before attaching) and 
x " — v{z n ) + d, and finally adding a new motor with 
its head at x n and its root at z n . Again, all the y- values 
have to be recalculated as stated by Eq. (3.1). Expressing 
y(z n ) through the positions of the neighbors (the index 
"— " describes the first attached motor on the left and 
"+" on the right hand side) yields 

-(z+ - z n ) + y+{z n - zJ) 



y 



d = 



(3.2) 



B. Equivalence to the model with a homogeneous 
force 

In the model described by now, the external force acts 
on one end of the backbone. This leads to some difficul- 
ties, e.g. one can only consider a semi-infinite chain with 
one boundary condition. Furthermore, the resulting so- 
lutions are not translationally invariant since the strain 
decreases along the backbone. 

Replacing the point force by a homogeneous one acting 
on the whole backbone with a density / = dF/dz would 
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allow us to perform the calculation on an infinite chain 
of completely equivalent motors. The ability to use peri- 
odic boundary conditions in the Monte-Carlo-simulation 
would be an additional advantage. 

Fortunately, both models, i.e. with a point force and 
a homogeneous force are actually equivalent. It is in- 
structive to show this equivalence in the continuum for- 
mulation of the model. Instead of the discrete set of 
variables yf we use a function y p {z). The continuum 
representation of Eqn. ([0])) is given by the following set 
of equations 



dz 2 



= -^5{z- z t )k(x[ - y p {z l )) 



7 



dyP_ 

dz 



F . 



(3.3a) 



(3.3b) 



z=0 



The first equation expresses the constant tension between 
the attached motors with jumps at the positions where 
the motor roots are placed. The second equation de- 
scribes the strain at the boundary of the backbone. 

On the other hand, in the homogeneous force model 
the strain grows linearly with z 



d 2 y 



dz 2 



f-Y,S(z-z l )k(x^y H (z l )) 



(3.4) 



H 



in 



f(NL- Zl ) 2 



(3.5) 



with F — NLf preserves the properties of the model. 
The transformation is shown schematically in Fig. |t]. Af- 
ter having shown the equivalence of both models, we 
can return to t he orig inal formulation. The transformed 
equations ( |3.l| ) and (3.2) become (in the following we 
omit the index H ) 

k{%i ~ Hi) - 1- — : \-Jz T ~ / = = 



Zi-l 



Zi+l — Zi 



(3.6) 



and 



V-{z+ - z n )+y+(z n - z-) 



+d-^-(z+-z n )(z n -Z-) 



(3.7) 



The additional term represents the displacement of a uni- 
formly loaded spring, tightly bound at its ends at z_ and 
z+. Of course, these equations also follow directly from 

(£j. 



C. Effective compliance of a semi-infinite chain 
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FIG. 7. Transformation from the model with a point forc e 
to the model with a continuous force, given through Eq. (3.5), 



on a typical configuration during the motion. The circles rep- 
resent attached motors, the detached ones lie on the line, z 
denotes the position a motor would have on an unstrained 
backbone, y its root position relative to z (Fig. ^|). 

Since we are dealing with a force constant in time and 
with quasistationary solutions, we can show that both 
models are equivalent up to a transformation which shifts 
the heads and roots of motors, depending on their z posi- 
tion. As may be easily verified by comparing Eqns. ( |3.3| ) 
and (3.4), the following transformation 




FIG. 8. Effective spring constant of a semi-infinite chain 
with randomly distributed displacements between bridges. 

A quantity frequently needed during the analytical so- 
lution of the model described in Sect. || is the elasticity 
of a semi-infinite stochastic chain as shown in Figure ||. 
It is defined as 



1 



dyi 
dF 



(3.8) 



where y -\ is a part of the solution of Eq. (3.1 
in Sect 
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As defined 

the spring constant of the motors is given by 
k. The values of U are distributed ra ndo mly with average 
L and the exponential distribution (|2.3| ). 

The compliance of a chain with a given configuration 
(given /^-values) can be calculated recursively. The chain 
is built up of a spring with stiffness k, connected in paral- 
lel to two other springs, which themselves are connected 
in series. The first one describes the piece of backbone 
with elasticity 7/Z1. The second spring is again a re- 
placement for another semi-infinite chain starting with 
the second motor. We denote its stiffness as K' 



G 



K e = k + 



h j 1_ 

7 ' A'£ 



(3.9) 



Repeating the same procedure for K' e etc. and finally av- 
eraging over all configurations I2, ■ ■ ■ ) with their sta- 
tistical weights yields 



1 



dl lP {h) Jdl 2P (l 2 ) 



-1 k+- 



7 fc+ J- 



(3.10) 

The convenient way to solve this high-dimensional inte- 
gral, however, is by using the Monte-Carlo method. Its 
result is shown in Fig. ^. 

One possibility to give an analytic approximation for 
(i-TjT 1 ) is the use of a mean- field like theory by assuming 
that all displacements U are exactly equal to their mean 
value L. Then the above expression becomes self-similar 
and we obtain the recursion relation 



1 



1 



(3.11) 



to the fact that the long-range coupling makes the de- 
tailed distribution of U irrelevant. The case of low 7 
is trivial since the motors get decoupled and the stiff- 
ness of the chain is determined solely by the first motor 



D. Linearity of the force velocity relation 



For rigidly coupled motors it was already shown in ||] 
that the velocity is linear in the applied force / as well 
as in the step size d. In the following we present a simple 
argument why this remains valid for elastically coupled 
motors. We denote the number of currently attached 
motors by n and their positions on the backbone by a 
vector x containing the comp onen ts x\ . . . x n . From the 
structure of Eqns. ( |3.6| ) and (3.7) it is evident that the 
process of attachment can be described by the following 
equation 



c(* +1 ) = Ax« + 



Ul 



fL 



u 2 d 



(3.13) 



with the solution 



jL\=±a +4 jL_i 

K e / 2 7 VV kL 



(3.12) 



where the (n + 1, 7i)-matrix A and the (n + l)-vectors Ui 
and u 2 depend in a complex way on {zi} and j/Lk. The 
detachment of one head is described by another (n—1, n)- 
matrix A as 



X (i+D = Ax C» 




FIG. 9. The effective compliance of a semi-infinite chain 
k{l/K e ). Thick line: Mont e-Carlo result. Thin line: approx- 
imation given by Eq. (3. IS). Dashed line: the ratio between 
the Monte-Carlo result and the approximation. 



In Fig. H the approximation, Eq. ( 3.12| ) is p lotted 
against the (exact) Monte-Carlo solution of ( 3.10| ). The 
deviations lie below 10% over the entire parameter range. 
The approximation becomes exact for both very soft and 
very rigid backbones. In the case of large 7 this is due 



(3.14) 



After a series of consecutive attachments and detach- 
ments this gives 



(3.15) 



Finally we set x^ ^ = and calculate the mean motor 
position (x) = trx/n. Since the time needed for i steps 
is proportional to t a , we obtain the relation 



ta V k 



(3.16) 



which is linear in d and in /. The constants C\ and Ci 
get independent of the mean number of motors for 
N -> 00. 

The force- velocity-relations for the elastic and for the 
stiff backbone are compared in Fig. ||. Due to the lin- 
earity of the force-velocity relation the problem can be 
separated in two parts: determining the backbone veloc- 
ity without external forces (/ = 0) and the velocity with 
external force but without power strokes (d = 0). The 
remaining work consists of determining the constants C\ 
and Ci. 
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E. Master-equation 

An adequate description of the model is given by 
P(. . . , I2, h; x, t; ri, r2, ■ ■ ■), the probability density to 
find a head at position x and with distances l\ = Zi~ Zi-i 
and n = Zi+i — Zi to its nearest attached neighbors, the 
distances I2 = 2^-1 — Zi^2 and r 2 = ^+2 — Zi+i between 
the nearest and the next nearest attached neighbors etc.. 
Of course, this distribution varies with time. Because the 
problem is linear in x and y, there is no need for deter- 
mining the correlations between the positions of different 
motor heads on the track. P can be expressed by 

P{. . . ,l 2 ,h;x,t;r 1 ,r 2 , ■ • •) = 

P(x,t)„„i 2 ,i i;ru r 2 ,...p(- ■ ■ ,h,h;n,r 2 , ■ ■ •) ■ (3.17) 

The second factor describes the probability for a motor 
to have the distances l\, r±, I2, r 2 , . . . to its neighbors. In 
a steady solution it is given by Eq. (2.3) 



e — ii/L g— ri/£ e — i 2 /i 



L 



L 



(3.18) 



The first factor gives the distribution of the head posi- 
tions x for the given set of distances. Since the transition 
rates are constant (strain-independent), the first factor 
does not influence the second one. 

We describe the temporal development of P in terms 
of a Master equation. The detachment of motors is de- 
scribed as drain, attachment as source. Additional terms 
result from the fact that the attachment/detachment of 
a motor also changes the I and r- values in its neighbor- 
hood. 



y\2 



y,-2 





1'2 



x 12 X U 



FIG. 10. Attachment of a new head. 



X,2 



The detachment rate is equal to the probability density divided through the mean life time of the attached state 



r d (. . . ,l2,k;x,t;r 1 ,r 2 , ■ ■ ■) = ^~P{- ■ ■ > h, h; x,t; r x , r 2 , . . .) 



(3.19) 



Once the distribution of z-values ( p(. . . , l 2 , h', 7"i, r 2 , ■ ■ •)) is m equilibrium, the attachment and detachment rate 
integrated over x have to be equal. They can only differ in the x-dependence. Thus we write the attachment rate the 
same way as the detachment rate except for a different factor containing the distribution of ^-positions of the newly 
attached heads. 



r a (. . . ,l2,k;x,t;r 1 ,r 2 , ■ ■ ■) = ^-p(- ■ ■ , h, h) ri,r 2 , . . •)- p "0M)... ) Wi;ri,r 2 , 



(3.20) 



The x-distribution of attaching motors depends on the ^-distributions of all the neighbors. It is determined as the 
integral over all (properly weighted) configurations which lead to a motor attachment at position x n : 



P n (x,t). 



dXnP(xn,t)... t l 3t l 2; l liri ... / dx rl P(x r i,t),,, i i uri . r 



.S(x n (. . . ,xn,x rl , ...) - x) 

(3.21) 



As follows form Eq. (3.7), x n can be expressed through 
yn and y rl 



xn = mn+ynii d _f_ lr 
h + n 27 



(3.22) 



which again a re fu nctions of all x% and can be determined 
through Eq. (jO]). The full Master equation is given in 
Appendix 



F. Zero load backbone velocity 

In the special case of zero external load (/ = 0) one can 
see that as long as the expectation value (2c}...,l a ,ii;ri,r 2 ,... 
is independent of the distances k and r i: i.e. 



JdxxP(x, 



t). 



,fa,h;rt,r 2 , 



= (x) , 



(3.23) 



the same holds for y^ (de term ined from ( |3.6| )) and for x n , 
which follows from Eq. ( |3.22 ) 



(x n ) = (y)+d= (x) + d 



(3.24) 



In other words, if the average head position of the exist- 
ing attached motors is not correlated to the distances be- 
tween them, the position where the head of a new motor 
attaches is uncorrelated too. Thus we have shown self- 
consistently that (x), ,.,i 2 ,i, :n ,r?,.. . do es not depend on U 
and Ti. The equation of motion (A2) for the expectation 
value of x simplifies to 



8 



« = 4<s> = ~((s n >-(s>) = r 



dt 



This means 



Ci 



1 



(3.25) 



(3.26) 



for the first coefficient in Eq. ( B . f t ) . The interesting point 
in this result is that it is independent of the backbone 
elasticity. As long as there is no external force, the ve- 
locity remains the same as in the case of a rigid backbone, 
this is v(f = 0) = d/t a . 



G. Slope of the force-velocity curve (correlations 
neglected) 

In the previous section we have shown that in the case 
of zero external force the average ^-position of a motor 
stays uncorrelated to the distances to its attached neigh- 
bors, which made the calculation of the force- velocity 
relation pretty easy. In the case of nonvanishing exter- 
nal force s (behavior described by the coefficient C2 in 
Eq. 3.16), the correlation doesn't vanish any more. How- 



ever, as a first approximation we may still try to neglect 
it. Later we will take correlations with distances to the 
nearest neighbors into account and show that they im- 
prove the result significantly. 

The solution is analog to the previous section, except 
that we set d = and / ^ in Eq. Q3.22Q . Another 
difference is that the average y- value differs from the av- 
erage a;- value by the mean force per attached motor (fL), 
divided through the motor stiffness k 



(V) 



(3.27) 



Instead of this rather intuitive argume nt t his equation 
can also be derived directly from Eq. (3.6) by summa- 
tion over i. 

This equation describes the y-position, average d ove r 
all motors. However, the quantities needed in Eq. ( 3.22j ) 
are the expectation values of y r \ and yn (motors at the 
edge of a gap whose distance to one neighbor is l\ + r*i 
and to the other neighbor randomly distributed). Even 
with uncorrelated x- values the expectation values (yn) 
and (2/7-1) are not the same as (y). 

For a gap width equal to the mean distance between 
attached motors both quantities will not differ. Other- 
wise the average value of the y-positions beneath a gap 
of width l\ + n differ from (y) by the excess force acting 
on the gap between them (f(l\ + r i)) compared to the 
force on the average gap (fL), multiplied by the effective 
compliance of a semi-infinite chain (Sect. Ill C ) 



(yn) = (yn) = (y) 



h + n-L 



(3.28) 



The denominator 2 describes the fact that the force is 
distributed equally to both edges. 

The mean x-position of the n e wly a ttache d mot or fol- 
lows by combining Eqns. (|3 22j) , (ft .271 ), and CH) : 

1 n\ 1 \ t ( I 1 \ h+ri- L L , lm ' 

{x )l ^ = {x) - f {\K;) — 2 — + k 



2l 



(3.29) 



Since we want to neglect correlations between (x) and 
(li, ri), we average over l\ and r\ ((l\) = (r\) = L) 



(3.30) 



In analogy to the previous section the final result for the 
dimensionless coefficient C2 is 



C 2 



k I I 



2\K, 



1 



kL 
2^ 



(3.31) 



The result is shown in Figure ^|. In the limit of very 
soft backbones (7 <C kL) the stall force becomes /stall = 
2d r y/L 2 , which is independent of the motor stiffness k. 
Our Monte-Carlo simulations show this behavior, how- 
ever with a different prefactor, which results from the 
correlations which were neglected in this approximation, 

/stall = —ft (3.32) 

with v « 1.64. 



H. First-order correlations 

Contrary to the previous subsection where we ne- 
glected the correlation between the position a; of a motor 
and the distances to its neighbors by using the ansatz 
([ 3.23 ), we now extend the calculation by taking corre- 



lation with distances to nearest neigh bors into account. 
We replace the approximation (3.23) by introducing a 
function describing these correlations 

dxxP(x,t)... t i 2t i i;ri , r2t ... = (x(t)) - ^/i (^j-; j- 

(3.33) 

The function /1 is scale invariant. It has to fulfill the 
condition 



dh I dr tl >Ur.n) t i[ |;^-) = 0. (3.S4) 



This also means that the head positions of the motors 
limiting a gap of width l\ + n already differ from (x) 



fL 

(x) h+ri , = (x) - — I dr 2 p(r 2 )v 



h + n r-2 
L 5 L 



(3.35) 



The equation ( 3.28 ) has to be extended by a term de- 
scribing this influence. Since the roots of the motors are 
connected to their heads via elastic constants k and to 
the rest of the track via effective constants K e — k, the 
correction in y corresponds to the x-correction attenu- 
ated by the factor k/K e . Finally, the refined Eq. ( 3.29] ) 
reads 



h + n-L 



L 



2 7 



(3.36) 



The correlation to farther neighbors is stil l neglected. This is expressed in the simplified equation of motion which is 
obtained by integrating both sides of (A2) over all Xi> 2 and pi> 2 



4 (x(t)} p(Xi; pi) = ^~ 
dt t a 



dh 



dr-i 



■ p(- ■ -,h;ri, . . .) 



(3.37) 



- ((x(t)) - (fL/k)f, (h/L; n/L)) (S(X 1 - h)6( Pl - n)) 

+ ((x(t)) - (fL/k)» (l 2 /L; h/L))(~5(X 1 - l 2 )S( Pl - h) + S(X 1 - l 2 )S( Pl - (h + n))) 
+ ((x(t)) - (fL/k)fi (r x /L; r 2 /L))(-S(X 1 - r l )5(p l - r 2 ) + 5(X 1 - (h + ri ))8( Pl - r 2 )) 
+ (^ n (*)) ii;ri (S(X 1 -h)S( Pl -r 1 )) 

+ ((x(t)) - (fL/k)p, (h/L; (h + n)/L) (-*(Ai - h)S( P i - (h + n)) + S(Xi - h)S( P i - h)) 
+ ((x(t)) - (fL/k)n ((h + n)/L; r 2 /L))(-S(X 1 - (h + n))6(pi - r 2 ) + S(X 1 - n )S( Pl - r 2 ))] 



which leads to an integral equation for p and v. Its scale invariant form, using the coefficient C 2 — vt a k/f, is 

A 

/ / 1 \ f\ 4- n - 1 \ 

1 

\ \K e / \ 2 ) 





C 2 = -3p(X:p) + J p(X;a)da + J p(a; p)da + ^fc y^-) (~ — ^ — ~ 



Xp kL 
~2~ 



+k 



1 



a p(X + p; a)da — (X + p)p(X; p) + / e a p(X; p + a) da + / e a p(X + a; p)d> 



(3.38) 



A self-consistent solution that holds for all A and p can be calculated numerically. The r esulti ng coefficient C 2 is shown 
in Fig. H and a typical shape of the /x-Function in Fig. 11. The correction v, defined in (3.32), gets the value v « 1.50. 
Taking first-order correlations into account improved the agreement between theory and simulation significantly. The 
remaining deviation, due to neglected correlations with further neighbors, is about 9% at low 7 and lies below 14% 
over the entire parameter range. 




FIG. 11. The function p(X; p) at 7/feL = 1. 



IV. STRAIN-DEPENDENT DETACHMENT 
RATES 



A. Stiff backbone 

A backbone with infinite stiffness (7 — > 00) means that 
the motor roots always keep their relative positions to 
each other, j/j = y for all motors i. In this case the posi- 
tions of m otor roots Zj naturally play no role. The Master 
equation (Al) for the z-independent probability density 
P(x, t) reduces to 



dt P ^ 



-r d (x) + r a (x) 



(4.1) 



with r d (x) = P(x)/t a ( x~y) rJ x) = S(x - y - d)L/Lt d , 
as follows from Eqns. ( 3.19^^22 ). The backbone position 
y is determined from Eq. (|3.6|), summed over i: 



k (x- y)P(x)dx = fL 



(4.2) 



Note that the norm of the distribution P gives the mean 
number of attached motors per length L, which, contrary 
to previous sections, i s no t necessarily equal 1. Quasis- 
tationary solutions of ([4.l|) are found with the ansatz 



P(X, t) = $(Z - y) 



y = vt, 



leading to 



(4.3) 



(4.4) 
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This equation is analytically solvable. For v > the 
non-divergent solution is 



6(d-0- (4.5) 



Finally, the force-velocity relation can be obtained 
through (Q f = kj £$(£)d£/L. The stall force is given 
as 



f( v = 0) = fed 



(4.6) 



For simplicity reasons, we model the strain-dependence 
of the detachment rate as an exponential function 



t a (0 = t a exp(a£) 



(4.7) 



The form of the forcc-velocity-relation, which was first 
measured by Hill Q, is fitted perfectly with ad = 0.58. 
Therefore, we use ad = 0.5 in the following. 



B. Elastic backbone 



In Section III we showed that for very flexible back- 



bones (j/kL <C 1) and strain-independent reaction rates 



(a = 0) the zero- load velocity remains 4-, while the 
stall force is limited through the backbone stiffness as 
f(v = 0) = 2jd/vL 2 . The mean strain on a motor is 
given as (£) = fL/k, or 2~/d/iskL at maximum load. For 
very flexible backbones this means a (£) -C 2ad w 1. 
The strain-dependence becomes negligible and the results 
from Sec. Ill are exact. From this simple argument we ex- 
pect that the force-velocity relation becomes linear. The 
crossover from the hyperbolic to the linear shape takes 
place at j/kL « 1. While the stall force (force at zero 
velocity) /stall = e ad kd/L is limited through the motor 
stiffness k for stiff backbones, it depends solely on the 
backbone stiffness for soft backbones, /stall = 2^d/vL 2 . 
This behavior is in agreement with our Monte-Carlo sim- 
ulations, the results are shown in Fig. ^. 
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APPENDIX A: MASTER-EQUATION 



In Eq. |3.17 we separated the probability density into a part depending on the distances between the z-values (U 
and Ti) and another part containing the x-positions of motors. Then we determined the attachment and detachment 
rates ( |3.19 3.20| ). Now we describe the temporal development of P in terms of Master equations. The detach- 
ment/attachment of one specific site leads to destruction respectively creation of the following states (see also Fig. 

My- 



Detachment 

destroyed 


Attachment 
created destroyed 


created 


(. . .,l 2 ,h;x;n,r 2 , ■ ■ ■) 
(. . .,h,ri;x r i;r 2 ,r 3 , . . .) 


(. ..,h,h;xn;h + ri,r2,-. ■) (• ■ ■ , h, h] xn; h + n,r 2 ,.. .) 
(. . . ,l 2 ,h + n.; avi ;r 2 ,r 3 , . . .) (■■■,l 2 ,h + n;2Vi; r 2 ,r 3 , ■ ■ ■) 


(. . .,l 2 ,h;x n ;n,r 2 , . . .) 
(. . . ,h,l 2 ;xn;h,n, . . .) 
(. . . ,h,ri;x r i;r 2 ,r3, . . .) 



The transition rates are given by Eqns. ( 3. IE , 3.2C ). The Master equation for P is 



dx I dli / dri 



t )...MM;Pi-P2-...P(- ■ ■ , Al! Pi: ■■■) z 
r d (. . .,l 2 ,h;x,t;ri,r 2 , . . .)(-£(£ - x) (• • ■ <5(Ai - h)8(p 1 - n) ■ ■ •) 

+P(x, t)..., hMuri ...{- ■ ■ ■ 5(X 2 - h)S(X 1 - l 2 )8( Pl - h) ••• + •• • 5(\ 2 - ; 3 )<5(A 1 - l 2 )S( Pl - {h + r 1 ))5(p 2 - r 2 ) ■ ■ •) 
+P{x, t)...j u r 1 -r a ,r a ..X 6 ( X 2 - h)6(\i - n)5( Pl - r 2 ) ■■■ + ■■ ■ 6{\ 2 - Z 2 )<S(Ai - (h + n))8{ P i - r 2 )S( P2 - r 3 ) • • •) 



+r a (. . . ,l 2 ,l 1 ;x,t;r 1 ,r 2 , ■ 

-\-P{x 7 ^)... J / 35 / 2 ;U+r , i,r2... (~ 
-|--P(x, ^)... ) i 2 ,ii+ri;r2,r'3... (~ 



• ■ ■ 5{Xi - l 2 )S( Pl - {h + r 1 ))5{p 2 - r 2 ) • 
■■■6(X 2 ~l 2 )6(X 1 -(h+r 1 ))6( Pl -r 2 )- 



.S(X 1 ~l 2 )S( Pl -l 1 )S(p 2 -r 1 )- 
■6(X 2 -l 1 )6(X 1 -r 1 )6( Pl -r 2 )- 



(Al) 

Rather than in the distribution itself we are interested in the expectation value {x(t)) ^ j r ^ . Its equation of 
motion follows directly form the Master equation. 

^.( x ( t ))-MM;pi,P2,-P(- ■ ■ , x i',Pi, ■ ■ ■) = j- J dh J dn ■ ■ ■?(■ ■ ■ ,k,h;r 1 ,r 2 , . . .) 
- (x(t))...,i 2 ,h;r 1 ,r 2 ,...S{(; - x) (• • • 5{Xi - h)6( Pl - n) ■ ■ •) 

+ (x(t))..., h ,i 2 - k ,rt..X- ■ ■ ■ <*(A 2 - h)S{X 1 - l 2 )S( Pl - h) ••• + •• • S(X 2 - l 3 )5(Xx - l 2 )S( Pl - (h + r x ))8{p 2 - r 2 ) ■ • •) 
+ (x(t))...,i uri ;r 2 ,r 3 ...( ^( A 2 - h)S(Xi - n)5( Pl - r 2 ) ■■■ + ■■ ■ S(X 2 - l 2 )S(Xi - (h + n))5( Pl - r 2 )6{ P2 - r 3 ) • ■ •) 



(x n (t)). M . rur2> ..s^ -!)(••• «(Ai - h)S( Pl - n) ■ ■ •) 

(x(t)},„,i 3 ,i 2 ;i 1+ri ,r 2 .,.( 8{Xi - h)S(pi - (h + r l ))8{p 2 - r 2 ) ■ 

{x{t)),„,l 2 ,h+rv,r 2 ,r 3 ...{ ^(^2 - h)5{Xi - (h + n))S( Pl - r 2 ) ■ 



■S(X 1 -l 2 )6(p 1 -h)6(p 2 -r 1 )- 
■S(X 2 -h)6(X 1 -r 1 )8{ Pl -r 2 )- 



(A2) 
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This equation simplifies further if the distances U and 
n are distr ibuted according to their equilibrium distribu- 
tion (3.18), which is certainly the case after the motors 
have been running for some time. 

Usually one is looking for the quasistationary solution 
with 



dt 



( x (t)) ...,A 2 ,Ai;pi,p 2 .. 



does not de- 



In the special case when (i fl )...,j a ,Ji;r 1 ,r a , 
pend on Zj and r^, the equation simplifies to ( 3.25 ). Tak- 
ing first-order correlations i nto a ccount but neglecting 
the higher ones leads to Eq. ( |3 .37 ). 
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